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Abstract 

We calculate the continuum value of the deconfining temperature in units of the string 
tension for SU(4), SU(6) and SU(8) gauge theories, and we recalculate its value for SU(2) 
and SU(3). We find that the iV-dependence for 2 < iV < 8 is well fitted by Tjy/a = 
0.596(4) + 0.453(30)/iV 2 , showing a rapid convergence to the large- N limit. We confirm our 
earlier result that the phase transition is first order for iV > 3 and that it becomes stronger 
with increasing N. We also confirm that as N increases the finite volume corrections become 
rapidly smaller and the phase transition becomes visible on ever smaller volumes. We interpret 
the latter as being due to the fact that the tension of the domain wall that separates the 
confining and deconfining phases increases rapidly with N. We speculate on the connection 
to Eguchi-Kawai reduction and to the idea of a Master Field. 



1 Introduction 



In a recent letter we presented lattice calculations of the finite temperature deconfining 
transition in SU(4) and SU(6) gauge theories. These calculations showed that the phase tran- 
sition is first order and becomes more strongly so as N grows. Our continuum extrapolations, 
taken together with older SU(2) and SU(3) results, showed that the 0{1/N 2 ) corrections to 
the N = oo limit appear to be small, even for SU(2). Finally, our finite volume studies showed 
that the finite volume corrections decrease as iV grows and we speculated on the simple picture 
that this suggested for SU(iV = oo). 

In this paper we will provide more accurate and more extensive calculations that confirm 
and make more precise our earlier study. For SU(4) and SU(6) we extend the range of 
lattice spacings at which we perform the calculations, thus allowing more reliable and precise 
continuum extrapolations. In addition we perform some calculations in SU(8) to strengthen 
our control of the N —>■ oo limit. We also perform some SU(2) and SU(3) calculations, 
both to check that our methods produce results in agreement with older (but still accurate) 
calculations and to investigate the crucial volume dependence in more detail. In addition we 
perform dedicated string tension calculations for all the above SU(iV) groups, generally at 
a value of the lattice coupling very close to the phase transition. This again improves the 
accuracy of our calculation and, as a by-product, leads to significant changes to the older 
continuum extrapolations of the SU(2) and SU(3) deconfining temperatures. 

The plan of the paper is as follows. In the next section we briefly summarise the details of 
the lattice setup. We then describe how we will identify the location and nature of the phase 
transition. We follow this with the details of our lattice results. We focus on T c / y/cr, the value 
of the deconfining temperature in units of the string tension, on the order of the transition, on 
the (Euclidean) latent heat, where first order, and on the iV-dependence of the finite volume 
corrections. We finish with some conjectures about the nature of physics at N = oo, and with 
some conclusions. 

We have, in addition, studied how topological fluctuations and (fc-)string masses vary 
across the phase transition, how the bulk transition and the thermodynamic latent heat vary 
with N, and other aspects of the transition. These, together with more details about the 
technical aspects of our calculations, will be presented in a later publication [2|. For some 
comments on the interest of SU(iV — > oo) gauge theories j^j, some references to relevant lattice 
calculations, and some speculations [4 on the nature of the deconfining transition, we refer 
the reader to our earlier paper pQ. 

2 Lattice setup 

Our lattice calculations are entirely conventional. We use hypercubic, periodic lattices with 
SU(iV) matrices, Ui, assigned to the links /. We use the usual plaquette action 
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where U p is the ordered product of the SU(iV) matrices around the boundary of the plaquette 
p, and which appears in the Euclidean path integral as e~ s . This lattice action becomes the 
usual Yang-Mills action with 

2N 

P = — (2) 
9 

in the continuum limit. The coupling g 2 is the bare coupling and, on a lattice of spacing a, 
may be regarded as a running coupling, g 2 (a), defined at the length scale a at which the action 
is defined. So to decrease a we decrease g 2 (a) and hence increase (3 (once a is small enough). 
As we vary iV we expect [3J, that we will need to keep constant the 't Hooft coupling, A, and 
its inverse, 7, 

\{a) = g 2 (a)N ; 7 = ± = JL (3) 

for a smooth large- N limit. Non-perturbative calculations jSlE] support this expectation. 

A periodic lattice whose size is L 3 L t with L ^> L t may be used to calculate thermal 
averages of the field theory at a temperature 

a((3)T = i. (4) 

Lit 

We can therefore study the deconfining transition by fixing L t and varying (3 so that T passes 
through T c . For simplicity we shall refer to the 3+1 dimensional system as being at the 
temperature T given by eqnfllj), although this can be deceptive (see below) and one needs to 
remember that the D=3+l system is no more than a theoretical framework for calculating 
thermal averages in a given three dimensional spatial slice at this value of T. 

The simulations are performed with a combination of heat bath and over-relaxation up- 
dates, as described in E]. We update all the N(N - l)/2 SU(2) subgroups of the SU(iV) 
link matrices. 



3 The phase transition 

We begin with some general comments about the deconfining phase transition and then we 
discuss how this phase transition appears in a finite volume and how one can extrapolate to 
the desired V — 00 thermodynamic limit. The discussion here will be heuristic, but can be 
formalised through standard finite size scaling theory, both for first [7J and for second |BJ order 
transitions. (We will summarise relevant aspects of this approach in our forthcoming paper 

3.1 General remarks 

The deconfining phase transition turns out to be second order for SU(2) and first order for 
SU(iV > 3). For a second order transition we have a diverging correlation length as T — > 
T c while for a first order transition various quantities will suffer a discontinuity at T = T c 
(including the correlation length which in general will be finite). This clear-cut distinction 
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only holds in the V — oo thermodynamic limit. On a finite volume, as in our numerical 
calculations, the two types of transition can be confounded and so a careful finite volume 
study is obligatory. We shall describe later on in this Section how we analyse the two kinds 
of phase transition. 

We are of course interested in the continuum limit of this phase transition. The continuum 
limit of the lattice theory is described by a 2'nd order phase transition at f3 = oo. This means 
that when we study the first order deconfining transition on the lattice, it will always appear to 
be weak; for example the latent heat is 0(a 4 ) and therefore small. Even on the coarser lattice 
spacings that are close to the bulk transition, which is first order for larger N, its latent heat is 
typically 1/50'th of the latter's. This is because the deconfining phase transition is a physical 
transition that only affects fluctuations on physical length scales, and leaves untouched the 
short distance fluctuations that usually dominate the values of composite operators. A typical 
lattice transition, like the bulk transition, affects fluctuations on all length scales. So to judge 
whether the first order phase transition is strong or weak we should compare the latent heat 
to a physical quantity, such as the T = gluon condensate. 

The phase transition is characterised by the fact that for T < T c the confining string 
tension is non-zero, a ^ 0, while for T > T c we have a = 0. To calculate the string tension 
we calculate the free energy of two static sources in the fundamental representation placed 
a distance r apart, F S g(r). Inserting such minimally coupled static sources into the action is 
equivalent to calculating the correlator of two Polyakov loops (ll{r)l p (0)) where 

n t =L t 

l p {n s ) = [ [ U (n s ,n t ). (5) 

n t =l 

(Here we introduce the notation, U^n), for the matrix on the link emanating from site n in 
direction fi.) In the confining phase at temperature T = l/aL t 

e -E^ll = ^j( r = aTlg )l p (0)) r ^°° e ~aMT)n z L t T <T C (6) 

where, to be definite, we take correlations in the z-direction. Since the operators l p are 
localised in a given z-slice, this correlator can be expressed in terms of the eigenstates and 
eigenvalues of the transfer matrix T z that translates us in the z-direction: 

(fj(r = an z )l p (0)) = J2 l(^c|/»| 2 e-~ (7) 

n 

where T z \n) = e~ arrin \n). If the lightest state in the sum is mi then 

(Zj(r = an z )l p (0)) " e"™ (8) 

We see from eqns ()6l8|) that 

ami = a 2 a{T)L t (9) 

and the fact that it is oc L t indicates that this state of the z-transfer matrix is the lightest 
flux loop that winds once around the t-torus. Once we have determined that the interesting 
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quantity is this lightest flux loop, we can alter the operator in the correlator to make the 
calculation more efficient. For example, we smear/block E] the links to have a better 
overlap onto the extended lightest state and we take p x = p y = sums of these smeared 
Polyakov loops so as to exclude explicitly the more energetic states with p ^ 0. 

In the confining phase the Polyakov loop cannot have an overlap onto a non-winding state 
since otherwise we would lose the factor of L t in the exponent of eqn(|fi|). This constraint is 
enforced by a C SU(N) symmetry of the action and integration measure: 

U (n s ,n t ) -> zU (n s ,n t ) Wn s (10) 

where z G Z^ is any one of the iV'th roots of unity and n t is some fixed value that one 
can freely choose. This symmetry ensures that (l p ) = since l p — > zl p under the symmetry 
transformation. 

In the deconfined phase cr(T) = and we therefore lose the factor of n z in the exponent 
of eqn(jnj) and so, in general, lim r .^ oo (Zj(r)Z p (0)) ^ 0, i.e. (l p ) ^ 0. In terms of the Z N 
symmetry described above, this is only possible if that symmetry is spontaneously broken. 
This tells us that for T > T c there are N degenerate deconfined vacua characterised by 
(l p ) = ce l2kn / N ; k — 0, N — 1 for c some real positive number. At very high T, where by 
asymptotic freedom g 2 (T) <C 1, one can calculate the tensions of the domain walls separating 
these vacua in perturbation theory jHj. It is important to note that what we actually know 
is that this Z^ symmetry breaking occurs in the D=3+l lattice system as one passes from 
l t = aL t > 1/T C to l t < 1/T C , and indeed in the continuum limit of that system. It is not 
however at all clear that it is actually a feature of the field theory at high T. That is to say, 
no-one has (we believe) succeeded in showing that this symmetry breaking is encoded in any 
thermal averages in a way that displays the existence of such Z^ bubbles and domain walls 
in the hot field theory 10 . To avoid tedious circumlocution, we shall usually speak of our 
D=3+l system as being at a temperature T = l/aL t although it is actually a framework 
within which to calculate thermal averages in a given spatial slice, and so one must always be 
careful when ascribing the detailed dynamics of this 3+1 dimensional system to the hot field 
theory. 

In summary, the natural order parameter for the deconfining phase transition is the average 
Polyakov loop, (l p ), and the mass of the flux loop that winds around the i-torus is the natural 
mass to focus upon, since it is simply proportional to cx(T). 

3.2 Characterising first order transitions 

If the deconfining transition is first order, the finite-T effective potential will have both a 
confining vacuum and N degenerate deconfined vacua. At T = T c all of these vacua are 
degenerate. On a finite volume there will be a finite range of T around T c in which the 
system has a significant probability to be in both confining and deconfined vacua, occasionally 
tunnelling between them. As we move away from T = T c , or (3 = /3 C , the degeneracy lifts and 
the confining/deconfined vacuum energies will split by some AE. The relative probabilities 
of the two vacua are given by e~ l3AE . Since AE oc V we see that the finite volume transition 
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region is typically 5(3 ~ 0(1/V). To calculate the V — oo value of (3 C , and hence of T c , we 
must first make a choice of (3 C on a finite volume, so that we can extrapolate it. A sensible 
criterion will use a value of (3 C (V) that is in the transition region, and therefore we generically 
can expect /3 C (V) = (3 c (oo) + 0(1 /V). We shall see examples of this below. 

In an infinite volume the average plaquette will jump across the transition and the lo- 
cation of this jump defines for us the (pseudo-) critical coupling (3 C (V = oo) and the critical 
temperature a((3 c )T c = 1/L t . In a finite volume the change will be continuous but rapid (if V 
is not too small). A useful quantity to characterise this behaviour is the specific heat C((3): 

j 2 C{$) = ^(u p ) = N p (u 2 p ) - N p (u p ) 2 (11) 

where u p = ReTrU p /N and u p is the average value of u p for a given lattice field, and N p = 6L 3 L t 
is the total number of plaquettes. The maximum of the specific heat is the maximum of 
d(u p )/df3 (neglecting the weakly varying (3 2 factor) and so the (3 at which it occurs provides 
at least one sensible definition of (3 C (V). 

Let u PiC and u p ^ be the average plaquette value for a lattice field in the confined and 
deconfined phase respectively, at (3 — (3 C (V). This value will fluctuate around its average 
value, (u PtC ) or (u P) d), by an amount that is 0(1/ y/N p ). As V — > oo we can neglect these 
fluctuations and we can also neglect the very rare fields where the system is tunnelling from 
one phase to another. In that case it is easy to see that the maximum of C {(3) occurs when 
the lattice fields are equally likely to be in the confined and deconfined phases, and in that 
case 

^JL-Cfa) = \({u p , c ) - (u p4 )) 2 = \l\ (12) 

where Lh is the latent heat per site of the phase transition and (3 C denotes the value of (3 at 
which the maximum of C(f3) occurs. (Note that this latent heat is that of the transition of 
the D=3+l system; it is related to, but differs from, the latent heat of the finite-T transition 
of the field theory.) The V^-dependence in eqn ffT^j) 

lim C(J3 C ) oc V (13) 

V — >oo 

is a direct consequence of the discontinuity in u p and so is characteristic of a first order 
transition. 

At a finite but large value of V, u p will fluctuate around (u p ) by an amount 0(1/ 'y/N p ) = 
0(l/y/V), i.e. (u p ) — (u p ) 2 = 0(1/V). (We can neglect the exponentially suppressed tunnelling 
fields.) It is then easy to show that C((3^) will receive a correction that is 0(1/V), as long 
as the fluctuations in the confined and deconfined phases differ. Thus we should expect the 
leading large-volume correction to both C{(3^) and f3^{V) to be 0{1/V). 

As we pointed out in Section I3~H the natural order parameter for the deconfining transition 
is the average Polyakov loop, (l p ). On a finite volume however, (l p ) = VT because of the 
tunnelling between the N deconfined vacua in the deconfined phase. (This is simply the 
statement that one only has spontaneous symmetry breaking in an infinite volume.) To 
finesse this problem we follow usual practice and replace (l p ) by (\l p \) where l p is the average 
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of the Polyakov loop over a single lattice field. Clearly (\J p \) ^ VT > T c . The price one pays 
is that now (\l p \) 7^ also for T < T c . However it is 0(1/^1^) in the confining phase and for 
a large volume this is a small price to pay. We can use this to define a quantity analogous to 
C(f3), where we replace the plaquette with the Polyakov loop. This is the normalised Polyakov 
loop susceptibility xi 

y = (\l P \ 2 )-(U 2 - (14) 

Just as with the specific heat we can define f3 c {V) to be the value of f3 at which xi has its 
maximum. At finite V the two definitions of /3 C (V) will differ by 0(1/V) terms, but will 
coincide in the thermodynamic limit. (We shall refer to them as /3* and /3f when we wish 
to distinguish them.) Again we have the characteristic first order behaviour xi V an d 
corrections down by one power of V: 

v ^v = ^ + v + - (15 > 

Just as for the specific heat, the maximum of xi occurs where the system is equally likely 
to be in the confined and deconfined phases (for V — > 00). One might think that a more 
natural criterion would be one where the deconfined phase was N times more probable than 
the confined one (simply because of the relative multiplicity of the former) and that it would 
therefore have smaller 0(1/ V) corrections. However we have not been able to construct a 
useful measure with such a behaviour to test this possibility. 

As we remarked above, we expect the V = 00 discontinuity to be smeared over an 0(1/V) 
transition region in a finite volume, so that definitions of f3 c (V) that use, for example, the 
peak of C(/3) or xi wm have corrections that are 0(1/ V") . In physical units, we expect the 
correction in the continuum limit to be 

r c (oo) - Tjy) y^oo h' 

T c (oo) VTl 

Inserting T C (V) = l/a([3 c (V))L t and Taylor expanding a(/3 c (V)) around V = 00 this translates, 
on an L 3 L t lattice, to the conventional finite size extrapolation jTTj 

&(V) - /3 c (oo) V =°° (17) 

where 

b! 2N 2 ti 



±\na(fi = p e ) ^lno( 7o )' (18) 

Thus if the leading finite-K correction is to be finite as V — > 00 we should find h oc A^ 2 . In 
PP we found some evidence that h is constant implying that the physical finite size correction 
vanishes as hi oc 1/iV 2 when iV — ► 00. Note that since hi is a dimensionless physical quantity, 
it receives lattice corrections that are 0(a 2 ). The value of h on the other hand also depends 
on how a varies with (5, i.e. on the violation of asymptotic scaling. 
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3.3 Characterising second order transitions 

If the deconfining transition is second order, the string tension er(T) will vanish smoothly 
as T — > T~ . The finite-T effective potential will have a single minimum for T < T c with 
a quadratic term that vanishes at T = T c . Thus at this temperature there is a diverging 
correlation length. This will be the inverse mass of the flux loop that winds around the t- 
torus; the vanishing at T c of this mass follows from eqnfJUJ) and from the vanishing of the 
string tension. For T > T c the former minimum at the origin becomes a maximum and a new 
minimum smoothly develops at a non-zero value of the order parameter (l p ). For SU(iV) there 
would be N degenerate such minima (as discussed in Section HQ} ; however we know that only 
SU(2) is second order so we can specialise to that case, in which case we have two deconfined 
vacua characterised by (l p ) oc ±1. 

Because the effective potential becomes flat near the transition, there will be large fluc- 
tuations. These will however be smooth, in contrast to the large tunnelling fluctuations of 
a first order transition. As T increases above T c one will begin to have distinct tunnelling 
transitions, but these will be between the deconfined vacua. 

Just as for first order transitions, C((3) or xi diverge at T c but the V dependence of this 
divergence is different and can be used to characterise the transition. Consider for example 
the specific heat. We may write it as 

J 2 C{P) = ^(u p ) = N p (u 2 p ) - N p (u p ) 2 = N p ((u p - (u p ))(u p - (u p ))) (19) 

Since the correlation length diverges as T — > T c , we expect that the correlation function of 
the plaquettes will vary as r oc°° e~ m °^ r where mo(T) — > as T — > T c and V — ► oo. Thus 

C{(3 C {V)) V oT V< (20) 

where £ < 1 depends on the critical exponents of the transition through ( = j/ud, using the 
standard notation. A similar derivation can be carried out for the loop susceptibility, xi- All 
this, as well as the ^/-dependence of /3 C (V), can be systematically analysed within the context 
of finite size scaling theory [S]. 

3.4 Tunnelling between phases 

On an L Lt lattice with Lt <C L we can view as follows the tunnelling that occurs between 
phases that have nearly the same free energy. An intermediate step required for tunnelling is 
that the L 3 periodic volume be split into two roughly equal domains each in one of the two 
phases. Separating these will be two domain walls. Since these walls have a finite tension the 
most probable configuration will be the one that minimises the area, so they will extend right 
across the short t-torus and across two of the spatial tori. The probability Pw(T) of forming 
such a bubble with two spatial walls of size I 2 = {aL) 2 is basically 

energy walls ^ a \v^ o„3_ t2t * 

P W (T) oc e ^ oc e — %~ oc e~ 2a awL Lt (21) 
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where aw is the tension of the wall (i.e. energy per unit area). 

We see from the L- dependence in eqn f!21|) that the probability of fields with such large 
domain walls is exponentially suppressed as the volume grows. This means that the 'time' 
the system spends in the transition between two phases is very small compared to the time it 
remains in one or the other phase: the transitions are rapid. It also means that the transitions 
are increasingly rare as L grows. These rare, large tunnelling fluctuations clearly distinguish 
between first and second order transitions once we are on a large enough volume. (Such 
tunnelling is also of course a feature of spontaneous symmetry breaking in a finite volume.) 

This provides a qualitative description of the transition between confined and deconfined 
phases at T = T c and between the N degenerate deconfined vacua for T > T c . In the latter 
case, for T ^> T c where g 2 {T) <C 1, one can use perturbation theory to evaluate the wall 
tension. One finds [9 for the wall tension between the k = and fc'th deconfined phases, 
characterised by l p ~ 1 and e l2k7T ^ N repectively, 

Att 2 T 2 

a° w k = k(N - k)^—(l + 0(g 2 N)) (22) 

where A(T) = g 2 (T)N is the 't Hooft coupling that one keeps constant for a smooth large- N 
limit. 

As we observed in jTj the deconfining phase transition appears to become more strongly first 
order as iV f, and the finite-V^ corrections appear to become rapidly smaller, e.g. h' oc 1/N 2 
in eqn (jl6|) . Our more extensive calculations in this paper will confirm this. The tunnelling 
near the phase transition can be seen on smaller V as N f, and for a given V (in physical 
units) the tunnelling becomes rapidly rarer as iV f. We interpret this as being due to the 
confining/deconfining interface tension, a c ^, growing with N. A reasonable guess is that for 
T ~ T c the tunnelling between the k = and k = N/2 deconfined phases (say N is even), will 
effectively pass through a sheet of confined phase so that 

1 k=0,k=N/2 AT 2 / O o\ 

a c4 oc -a w oc N (23) 

using eqn (|22j) . From eqn(|21|) we see that this means that the tunnelling has a probability 

P Cjd {L,T) oce"^- (24) 

so that to see tunnelling equally easily we need to reduce the volume L oc 1/N as we increase 
N. This roughly describes what we shall see in our calculations. 

The assumptions of our above argument can, in fact, be weakened considerably. All we 
need is that a w ~ 0,k ~ N ^ 2 ~ 0(N 2 ) for T ~ T c , and then the fact that cr^7°' fc_7V//2 < 2a c d obtains 
for us eqn(|23|) as a lower bound on the iV-dependence of the domain wall tension. We remark 
that the growth of domain wall surface tensions as a power of iV appears to be a common 
phenomenon (see SectionEJ) so it is no surprise that it also occurs in the case of the deconfining 
phase transition. 
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4 Results 



4.1 Setting the scale 

To set the scale of the lattice spacing in physical units we use the T = string tension, o. We 
extract o from the mass of the lightest flux loop that winds around the spatial torus, which 
can be obtained from correlations of smeared Polyakov loops as described, for example, in [Sj. 
We generally use spatial tori of length I = aL > 3/y/a, since earlier calculations [3] suggest 
that finite volume corrections are then accurately given by just the leading universal string 
correction, i.e. 

arrip(L) = a 2 aL (25) 

assuming the confining string is a simple bosonic string as suggested by [T3] . 

We have performed string tension calculations for all the SU(iV) groups discussed in this 
paper. The calculations are taken from a study of improved gluonic lattice operators [Hj and 
we use a type of blocking/ smearing that is much more effective than in earlier calculations. 
This leads to significant changes in earlier estimates of the continuum value of T c /^Ja for both 
SU(2) and SU(3). 

Our calculated values of a^/a are listed in Tableland Tabled Once we determine the 
critical coupling (3 c (oa) for a given value of L t , we interpolate our values of a^Jo to that value 
of (3 and this provides our lattice estimate of T c 

aT r T r 1 1 . , 

a(p c ) = — — . (26) 



a^Jo L t a{j3 c )^/a ' L t T c 

Repeating the calculation for various L t , we can attempt to extrapolate the results to the 
continuum limit using a leading 0(a) lattice correction 

T c (a) T c {a = 0) , _ 2 



vM a ) ^a(a = 0) 
where c is some unknown quantity that is fitted. 



+ caV (27) 



4.2 Locating the phase transition 

To locate a first order transition we search for values of /3 at which there is tunnelling between 
the confined and deconfined phases. (See e.g. Fig.l and Fig.6 of [1].) This is easiest to do 
on smaller volumes where the tunnelling is sufficiently frequent that it can be seen in modest 
Monte Carlo runs. (Although if the volume is too small the transition will be completely 
washed out.) One can then use the /3 C (V) extracted from a smaller volume as a starting point 
for calculations at larger V. 

We typically carry out O(200K) sweeps at each (3 and have 2 to 5 useful values of (3. They 
will be close enough together that the distributions of u p have a strong overlap, and we then 
use standard reweighting techniques [12] to produce Xi{&) an d C{(3) as continuous functions 
of (3 and to obtain an accurate value for the position of the maximum, (3^(V) or (3*(V). An 
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example is shown in FigO We then repeat the calculation for various volumes, L 3 , at a fixed 
value of L t and perform a finite size study (see below) to obtain /3 c (oo), using eqn(fT7|). the 
latent heat, using eqn fT2|) . and so on. We then proceed to the lattice and continuum values 
of T c /y/a as described in Section I4~T1 

In practice we carry out the finite volume study at one value of L t for each value of N and 
then use the parameter h obtained from a fit of the form in eqn()17|) at the other values of L t 
where we perform calculations typically on just one volume. This is the procedure used in 
past SU(3) calculations, e.g. jTT], but is dangerous because the corrections to h in eqnffTTj) are 
not merely 0(a 2 ), as they would be for h! in eqn([TB|). but include the violations of asymptotic 
scaling in the d\aa/df3 factor. This is a weakness of this calculation, which can only be 
remedied by a study on a much larger scale than ours. 

For the case of SU(2), where the transition is second order, it turns out that the specific 
heat C{j3) has no visible peak as we vary (3 and its value does not grow with V. This anomalous 
behaviour arises because the plaquette has almost no overlap onto any correlation length that 
diverges. This interesting behaviour will be discussed in more detail in f2J. Here we merely 
note that one can instead focus upon the Polyakov loop susceptibility xu as is usually done, 
and that exhibits all the characteristics of a second order phase transition. 

4.3 Finite size study at L t = 5 

We perform our finite- 1/ studies at L t = 5 for all the values of N that we investigate. For a 
quick global view of what is going on we plot in Fig.2 the value of Inxi (extracted at its peak 
value) against In V for all our values of N. We see that for N > 3 the values fall approximately 
on a straight line of slope unity. This is just the slope appropriate to a first order transition 
(see eqns ()13Jlf)j) ) and is in contrast to the SU(2) line that is flatter (with a slope of ~ 0.63) 
indicating that the transition there must be second order. 

To see the corrections, a plot of the normalised susceptibility maxima x™ ax /^ versus 1/V is 
more useful, and this we show in FigHfor SU(2) and SU(3). We see that for SU(2) xf/V 
as V — » oo indicating it is second order. (Our largest SU(2) lattice is L = 64.) For SU(3) the 
behaviour at smaller V is similar to that of SU(2), indicating a growing correlation length, but 
eventually it turns upwards, showing that it is indeed dominated by a first-order discontinuity. 
So SU(3) is 'weakly' first order. For N > 4 (FigHJ) the picture is much less ambiguous: there 
is the expected correction that is linear in 1/V with a flattening only at small V where 
presumably the correlation length is being felt. As N j the values of V where the flattening 
occurs is at ever smaller V indicating that the strength of the first order transition is growing 
(in the sense that the discontinuity dominates over the contribution of the longest correlation 
length at ever smaller V). 

In Fig|5]we perform similar plots for the normalised specific heat maxima C maiX /V. We see 
a much nicer linear behaviour that is easy to extrapolate. The reason for this is that, unlike 
Xu C{(3) is not sensitive to the lightest correlation length (which is proportional to a(T)). 
Indeed we do not show the SU(2) curve, because we cannot identify a maximum - and if we 
did its value would be independent of V (in our range of V). 

We can use the maximum of xi to define /3 C (V) and from a fit to eqn ffTTj) we can extract 
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h, the coefficient of the 1/V correction. This is tabulated in Table El We recall that the 
usually quoted SU(3) value is h < 0.1 fllj. We find it impossible to extract a reliable value 
of h from our values of f3 c {V) at L t = 5 (although we have some calculations for L t = 4 that 
agree precisely with earlier calculations) and we put that down to the weakness of the SU(3) 
transition - or perhaps a large statistical fluctuation in our calculations. In any case, it is 
already clear from Table El that we do not have h oc iV 2 as one would need for a finite h' in 
the iV = oo limit. (See eans (jlh117H8j) .'l Rather it looks as though h! oc 1/N 2 . 

4.4 Latent heat 

We see from eqn(JT2J) that we can obtain L h , the latent heat per plaquette, from the V — oo 
extrapolation of C max /V. Doing so for our L t = 5 calculations, we obtain the value of L h for 
each N and we plot the results in FigEl 

The first thing we observe is that Lh grows with N. Moreover the growth is consistent 
with being due to a 0(1 /N 2 ) correction, as shown by the dashed line in the plot. It is clearly 
no surprise that the latent heat vanishes for 2 < N < 3, as exemplified by our fit, so that the 
N = 2 transition is second order. 

The lattice latent heat is 0(a 4 ). Can a variation of a with N account for the variation 
we see in FigEP If we express a in units of T c then the answer is of course no: by definition 
a = l/L t T c = 1/5T C is fixed for all N. If instead we use a to set the scale, then a 4 does indeed 
vary as we go from N = 3 to N = 6, but only by about ~ 30% which is a very small part of the 
observed change. Thus we conclude that in physical units the latent heat increases strongly 
as we go from N = 3 to larger N, indicating that the first order transition strengthens as iV f. 

Finally we remark that if the latent heat goes to a non-zero constant as iV — > oo, this 
means that it is proportional to the gluon condensate (with the 't Hooft running coupling) 
[IE! so that it is 0(N 2 ). 

4.5 T c /Vgt for all N 

We list our values of /3 C (V = oo) in Table |U and Table El The SU(2) values include older 
calculations 17j as do the SU(3) ones [TT1 IT2*]. The (interpolated) string tensions are new and 
more accurate than those used previously. For our SU(4) and SU(6) calculations at L t = 6, 8 
we use the h extracted at L t = 5. However we double the error so as to (hopefully) cover any 
change of h with L t . For SU(8) we also extract h for L t = 5, but because the error on it is 
already very large, we feel no need to increase it at larger Lt. 

We use eqn (f2T)j) to obtain the values of T c j ^/a shown, and then use eqn (|2*7]) to obtain the 
continuum extrapolation, as tabulated in Table El (In the SU(4) case where the best x 2 is 
very poor, we use an error that covers the extrapolated values obtained using any two ft values 
out of the three.) 

Finally we take our various continuum values of T c /yja and plot them against 1/N 2 in 
Fig. [7| We see that the values are all well described by a leading large- iV correction of the 
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expected |3j functional form: 



0.596(4) + 



0.453(30) 

W 2 



(28) 



The coefficient is modest, indicating that for this physics at least, all values of N are close to 
N = oo, despite the fact that SU(2) is second order while larger groups are first order. Only 
the latent heat appears to have large corrections in N. 

4.6 F-dependence as N f 

The values of h in Table E2 strengthen our earlier observation PQ that h is roughly constant 
indicating that the finite volume corrections disappear as oc 1/N 2 . 

It is also clear from e.g. FigJ3] and 0] that the first order character of the transition becomes 
clearer on ever smaller volumes as V f. We have interpreted this as being related to the domain 
wall surface tension growing with N, perhaps as in eqns ()22|23|) . To analyse this further we 
have roughly calculated the number of sweeps that the lattice spends in each of the confined 
and deconfined phases when it is very close to (3 C (V). We call this number the persistence 
time, t p . From the statistical errors on our T = string tension calculations, we infer that the 
efficiency of our Monte Carlo in decorrelating physical fluctuations does not vary much with 
N, so that comparing r p across N is physically meaningful. We plot r p for various volumes 
and various N in FigJHl against L 2 . We see that it is consistent with being r p oc exp —(N a L) 2 
with a unity or greater. In fact the L dependence is not really well determined: a plot against 
L 3 looks equally good. If we simply take the largest L where we have roughly equal r v for 
iV = 3,4,6 i.e. L = 32,20,12 respectively, then a trend to be approaching a function of 
N x L for larger N is certainly consistent. However this analysis has to be regarded as very 
preliminary at this stage, and only indicative of what one might eventually achieve j2]. 

5 Some conjectures about N = oo 
5.1 Reduced models 

On a given spatial volume, V = L 3 , at fixed L t and at values of j3 close to the corresponding 
phase transition, we have found the tunnelling to be rapidly suppressed as N f. Roughly 
speaking if one wants to keep the tunnelling rate constant one needs to decrease the volume as 
L oc 1/N a , where a appears to be at least unity. This is roughly consistent with our conjecture 
that the surface tension of the confining/deconfining domain wall varies as a c ^ oc N 2 , and in 
any case that this provides a lower bound on its variation with N. We also saw that the shift 
of aT c {V) from aT c (oo) appears to be 0(1/N 2 ), suggesting that as N — > oo one can locate 
the transition on volumes V — ► 0. 

As we pointed out in [T| this cannot really be the case, because once L < l/a(/3)T c ~ L t 
the spatial version of the center symmetry in eqn()10j) will become spontaneously broken. For 
example if we are considering L x < L t then the relevant symmetry is as in eqn()10|) but with 
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U x replacing U , n x replacing n t , and n s counting the points in the (y,z,t) 3-slice. This 
breaking occurs for the same reason that the symmetry is broken in the temporal direction 
for L t < l/a((3)T c ; after all the label of the Euclidean direction does not matter. (For finite 
N this breaking will itself be smoothened by finite volume effects; however as iV — > oo the 
observed loss of finite- V corrections suggests it will become precise, even on a finite volume.) 
For this reason we conjectured in pQ| that as iV — > oo the deconfining transition will occur 
at exactly the same /3 C for all volumes L > L t = l/a(/3)T c , i.e. the deconfining physics is 
independent of the volume down to a critical volume, and we speculated upon the relation 
of this to Witten's idea of an N = oo Master Field that is necessarily translation invariant 
[IE]. We also note the relation of this to recent conjectures [THj about how the Eguchi-Kawai 
reduction [20] translates to the continuum. 

Of course, as iV — > oo the transitions on this minimal V will presumably become rapidly 
rarer and the transition will become ever harder to locate. Here we speculate on a way of 
reducing V below this minimal value. This is based on the observation that if we impose 
an {x, t} twist for T > T c , it induces a corresponding Zjq domain wall separating two of 
the Z N deconfined phases. That is to say, it disorders the Z N symmetry breaking. Indeed 
this domain wall can be interpreted [21] as a 't Hooft dual vortex in its 'Higgs' phase. If 
L x becomes small there is no room for the two phases and the wall, and we expect that 
the corresponding symmetry breaking is suppressed. Thus we conjecture that on an L 3 L t 
volume at /3 C , if we reduce L < L t but introduce a (maximal) twist in {x, y}, {x, z} and {y, z} 
we suppress the symmetry breaking and can continue to reduce L, perhaps as 1/N, while 
maintaining the physics of the deconfining transition. As N — > oo we can therefore study T c 
on an infinitesimal volume: in a sense on a single spatial point. And perhaps other physical 
quantities can also be so studied. This reduction connects naturally to the twisted Eguchi 
Kawai reduced model [201 [22] • 

5.2 Master field(s) 

The factorisation of gauge invariant operators at large- N, and the fact that fluctuations vanish 
in leading order, led to an old speculation ^H] that in the N = oo limit there is only a single 
Master field in the path integral (which is necessarily translation invariant for gauge invariant 
quantities). It has long been realised that this picture is deceptively simple (e.g. [23]) • 

In the case of interest here, we have seen that the fluctuations near T = T c between the 
confined and deconfined phases, are suppressed exponentially in iV 2 (at least) on a volume 
that is fixed in physical units as N f. If one wishes to talk of a 'Master Field' one has to talk 
of at least 2: one confined and one deconfined. (For the D=3+l Euclidean system there are 
presumably iV — > oo deconfined Master Fields, separated by infinite barriers in that limit.) 

We note that as T deviates from T c the difference in free energy density between the 
confined and deconfined vacua will be 0(N 2 ) and for large enough bubbles will outweigh the 
suppression of the domain walls (if that is also 0(N 2 ) as we have suggested). However the 
system has to pass through small bubbles to reach these large bubbles, and the probability 
of these will continue to be dominated by the exponential suppression of the surface area 
of the bubble. That is to say, one can imagine a scenario at iV = oo in which there is no 
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transition between confining and deconfining phases at any value of T - a kind of infinite 
super-cooling/heating. 

This picture of different phases separated by infinite barriers at N = oo is also char- 
acteristic of the lattice bulk transition Ej and of the multiple vacua inferred from the 
^-dependence of the SU(iV) gauge theory |2H EE]- All this suggests that the gauge theory 
possesses an infinite number of 'Master Fields' that mutually decouple precisely at N = oo 
but whose simultaneous existence ensures a rich physics of that limiting theory. 

6 Conclusions 

In this paper we have obtained some properties of the deconfining transition in SU(iV) gauge 
theories. We have done so for enough values of N that we feel confident in our control of 
the N —>■ oo limit. At the same time, our dedicated string tension calculations allow us to 
improve significantly upon older continuum extrapolations for SU(2) and SU(3). 

We find the phase transition is first order for N > 3, and that the latent heat (of the 
D = 3 + 1 Euclidean system) grows with N, as shown in Fig® showing that the transition 
becomes more strongly first order as N |. Indeed it is clear from FigJU] that the SU(3) 
transition is quite weakly first-order, and that a naive interpolation of the latent heat as a 
function of N predicts that the transition ceases to be first order somewhere between N = 2 
and N = 3. 

We also find that as N grows the phase transition becomes clearly defined on ever smaller 
volumes (see e.g. FigEJEJ) - We argue that this is due to the interface tension between the 
confined and deconfined phases growing at least as rapidly as N 2 . At the same time we 
find that the finite- V corrections to T c appear to be 0(1/N 2 ) so that at iV = oo T c can 
be calculated on small volumes. We speculated in Section how one might be able to take 
the limit V — > in a way that is reminiscent of the twisted Eguchi-Kawai model. We also 
speculated on how this might fit in with the old idea that N = oo is dominated by a Master 
Field that is translation invariant for gauge invariant observables. We plan to test the idea in 
Section IB~Tl not only for T c but for other physical quantities as well. 

The values of T c / yja that we obtain in the continuum limit show a modest iV-dependence, 
displayed in Fig. This dependence can be well fitted with a leading 0(1/N 2 ) correction all 
the way down to SU(2), as in eqn ()28|) . Such modest corrections are in fact characteristic of 
the N dependence of many physical quantities • 

As we have stressed here and in there are many ways in which our calculations need to 
be improved. Perhaps their main weakness lies in the control of finite- V corrections. Because 
such studies are expensive we have followed the usual strategy of performing a finite- study 
only at the coarsest value of the lattice spacing a and then using that to estimate finite-V^ 
corrections at other values of a. That this is likely to be dangerous is indicated by the nature 
of the theoretical expression, eqn (|TH|) . Another weakness concerns our latent heat calculation 
which, being coupled to the finite-^ study, has again only been performed for the coarsest 
value of a. Here in fact we can do better, using a different technique, and in j2] we will do so. 
In that paper we will also present calculations of the correlation lengths as T — > T c , obtained 
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separately for the confined and deconfined phases. This will provide another measure of the 
strength of the transition. In addition we will present some results for the physical latent heat 
of the gauge theory (which is different from the 'Euclidean' latent heat that we discuss in this 
paper), for the way topological fluctuations change across T = T c , and we will present the 
results discussed in this paper in greater detail. 
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SU(2) 


SU(3) 


P 


L 




P 


L 




1.8800 


4 


0.773(16) 


5.6925 


8 


0.3970(19) 


2.1768 


8 


0.5149(77) 


5.6993 


8 


0.3933(16) 


2.2986 


10 


0.3667(18) 


5.7995 


10 


0.3143(14) 


2.3726 


12 


0.2879(10) 


5.8945 


12 


0.2607(11) 


2.4265 


16 


0.2388(9) 


6.0625 


16 


0.19466(73) 


2.5115 


20 


0.17680(76) 


6.3380 


24 


0.12946(61) 



Table 1: SU(2) and SU(3) string tensions calculated at the indicated values of (3 from confining 
strings of length aL that wind around the spatial torus. 
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SU(6) 


SU(8) 
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a^Jo 


P 


L 




P 


L 


a^Ja 


10.637 


10 


0.3254(11) 


24.500 


10 


0.3420(19) 


43.85 


8 


0.3646(25) 


10.789 


12 


0.27015(86) 


24.845 


12 


0.2801(13) 


44.00 


10 


0.3406(20) 


11.085 


16 


0.19871(85) 


25.452 


16 


0.20992(82) 


44.35 


10 


0.2991(20) 














44.85 


12 


0.2596(24) 



Table 2: SU(4), SU(6) and SU(8) string tensions calculated at the indicated values of j3 from 
confining strings of length aL that wind around the spatial torus. 





h 


L E x 2 /n df 


SU(4) 
SU(6) 
SU(8) 


0.090(17) 
0.112(19) 
0.063(70) 


[12,20] 0.49 
[8,14] 0.04 
[6,8] 



Table 3: The coefficient h of the 1/V correction to P£(V) for the groups shown, obtained 
from fits (best x 2 shown) to the range of volumes shown. 
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SU(2) 


SU(3) 


In 


Pc 




Lt 


Pc 




2 


1.8800(30) 


0.647(15) 








3 


2.1768(30) 


0.6474(111) 


4 


5.69236(15) 


0.6296(30) 


4 


2.2986(6) 


0.6818(35) 


5 


5.8000(5)* 


0.6370(29) 


5 


2.37136(54)* 


0.6918(29) 


6 


5.8941(12) 


0.6388(32) 


6 


2.4271(17)* 


0.6994(50) 


8 


6.0625(18) 


0.6421(31) 


8 


2.5090(6)* 


0.7008(34) 


12 


6.3385(55) 


0.6442(51) 



Table 4: Critical values of P, extrapolated to V — oo, for SU(2) and SU(3) for the values of 
L t shown; with the corresponding values of the deconfining temperature in units of the string 
tension. Starred values are new calculations of P c . 





SU(4) 


SU(6) 


SU(8) 




Pc 


T c /y/a 


Pc 


T c /y/a 


Pc 


T c /y/a 


5 
6 
8 


10.63709(72) 
10.7893(23) 
11.0848(23) 


0.6146(23) 
0.6171(26) 
0.6289(31) 


24.5139(24) 
24.8458(33) 
25.4712(62) 


0.5894(36) 
0.5952(32) 
0.6009(29) 


43.978(22) 
44.496(3) 
45.606(35) 


0.5814(68) 
0.5807(71) 
0.5964(154) 



Table 5: Critical values of p, extrapolated to V = oo, for SU(4), SU(6) and SU(8) for the 
values of L t shown; with the corresponding values of the deconfining temperature in units of 
the string tension. 





lim a ^ T c /y/a 


L t e 


X 2 /n df 


SU(2) 


0.7091(36) 


[3,8] 


0.28 


SU(3) 


0.6462(30) 


[4,12] 


0.05 


SU(4) 


0.634(12) 


[5,8] 


2.05 


SU(6) 


0.6078(52) 


[5,8] 


0.00 


SU(8) 


0.594(20) 


[5,8] 


0.53 


SU(oo) 


0.5960(41) 




0.29 



Table 6: Continuum extrapolations of T c j ' yja with the range of L t values used and the x 2 P er 
degree of freedom of the fit. Also shown is the result of the iV — > oo extrapolation of these 
values, as described in the text. 
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Figure 1: Reweighted susceptibility \i as a function of (3 for SU(3) on a L = 32 lattice at 
L t = 5. The location of the maximum, f3 c , and the value of the maximum, x™ aX ' are gi ven as 
well. 
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Figure 2: The value of lnxf 1 ^ plotted versus In V for all our values of N at L t = 5. The full 
lines are best fits for N = 3, 4, 6, 8 with a slope close to unity while the dashed line is the best 
fit for N = 2 with a slope of ~ 0.63. 
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Figure 3: The normalised susceptibility maxima x™ ax /^ plotted against 1/V for N = 2,3 at 
L t = 5. The dashed line is a best fit with the leading scaling behaviour extracted from FigJ2] 
and a leading 0{1/V) correction. 
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Figure 4: The normalised susceptibility maxima x" iax /^ plotted against 1/V for N = 4, 6 at 
L t = 5. The dashed lines are the best fits with 0(1 /V) and 0(1 /V 2 ) corrections. Note the 
different scale for SU(6) on the rr-axis. 
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Figure 5: The normalised specific heat maxima C max /V plotted against 1/V for N = 3,4,6 
at L t — 5. The straight lines are best fits with a leading 0(1 /V) correction. Also shown are 
the V — > oo limits. 
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Figure 6: The latent heat, L h , plotted versus N at L t = 5. The dashed line is a large- N 
extrapolation with a leading 0(1/N 2 ) correction. 
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Figure 7: The deconflning temperature, T c , in units of the string tension, a plotted versus 
1/N 2 . The dashed line is the large-iV extrapolation with a leading 0(1/N 2 ) correction. 
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Figure 8: The persistence time, r p , plotted versus L 2 for various N at L t = 5. The straight 
lines are fits of the form r p oc exp —c(N)L 2 . 
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